Double Higgs Production and Quadratic Divergence Cancellation in Little Higgs 

Models with T-Parity 



Claudio 0. Ditfl 

Department of Physics, Universidad Tecnica Federico Santa Maria, Valparaiso, Chile 

Rogerio Rosenfelcf] 

Instituto de Fisica Teorica - UNESP, Rua Pamplona, 145, 01405-900 Sao Paulo, SP, Brasil 

Alfonso ZerweklJI 

Department of Physics, Universidad Austral de Chile, Valdivia, Chile 
(Dated: July 2005) 

We analyze double Higgs boson production at the Large Hadron Collider in the context of Little 
Higgs models. In double Higgs production, the diagrams involved are directly related to those 
that cause the cancellation of the quadratic divergence of the Higgs self-energy, providing a robust 
prediction for this class of models. We find that in extensions of this model with the inclusion of 
a so-called T-parity, there is a significant enhancement in the cross sections as compared to the 
Standard Model. 



I. INTRODUCTION 

The presence of quadratic divergences in loop corrections to the scalar Higgs boson self-energy is responsible 
for the so-called hierarchy problem of the Standard Model (SM) ; namely, there is no natural way of having a 
"light" mass (i.e. ~ 10 2 GeV) for the Higgs given that loop corrections induce contributions to the mass of the 
order of the GUT scale -or in general the high energy scale above which new physics enters and the SM ceases 
to be an effective theory. In supersymmetric extensions of the SM this problem is absent since the bosons are 
protected from quadratic divergences by their relation to supersymmetric (fermion) partners^. The hierarchy 
problem is also absent in models in which the electroweak symmetry is dynamically broken, since the scalar 
particles are not fundamental but composite in those cases {2|. 

Recently a new kind of model was proposed which can solve the hierarchy problem of the scalar Higgs boson. 
In these models, called Little Higgs (LH) models 0, the Higgs boson is a pseudo-Goldstone boson and its mass 
is protected by a global symmetry and, unlike supersymmetry, quadratic divergence cancellations are due to 
contributions from new particles with the same spin. 

The phenomenology of these models has been discussed with respect to indirect effects on precision measure- 
ments Q and direct production of the new particles introduced Q . 

Since these early contributions, several variations in the LH framework have been proposed 6]. However, the 
cancellation of quadratic divergences is inherent to any LH model and this requires definite relations among 
certain couplings. Therefore, any process that involves exclusively these couplings is a robust prediction of the 
LH mechanism regardless of model variations. 

In this article we study a process that has ingrained in it the cancellation of quadratic divergences of top-quark 
loops, namely double Higgs production. In section II we review the model and the derivation of the masses and 
couplings of the relevant particles by appropriate diagonalization of mass matrices and show the cancellation 
of quadratic divergences directly in the broken symmetry phase. In section III we derive the amplitudes for 
double Higgs production at the Large Hadron Collider (LHC) and compute the cross sections. Our results are 
presented in section IV and we conclude in section V. 
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II. MASSES, COUPLINGS AND QUADRATIC DIVERGENCES IN THE LITTLEST HIGGS 

MODEL 



There are many variations of Little Higgs models today, which differ in the symmetry groups and represen- 
tations of the scalar multiplets, but they all have in common a mechanism of cancellation of the quadratic 
divergence for the mass of the lightest remaining scalars at one loop order. After the spontaneous breakdown of 
a global underlying symmetry at a scale Awf (supposedly not much higher than a few TeV to avoid fine tuning) , 
the model contains a large multiplet of pseudo-Goldstone bosons, which includes the SM Higgs doublet. While 
most members of the multiplet receive large masses (again, of a few TeV), the mass of the Higgs boson is 
protected from quadratic divergences at one loop, and therefore remains naturally smaller. The cancellation is 
related to the existence of an extra (heavier) top-like quark and its interactions with the scalar sector, feature 
which is common to all LH models. Consequently, a good test to distinguish a little Higgs from other cases 
should be based on a signal sensitive to this particular feature of divergence cancellation and rather insensitive 
to other features. The Higgs pair production at LHC is one of such signals, since it is based on exactly the 
same diagrams that enter the quadratic divergence cancellation (Fig. 1), except for the insertion of two gluons 
(Fig. 2 and 3). 

In order to work out the details, we make use of the Littlest Higgs model, which is a simple case but contains 
all the necessary ingredients. After spontaneous breakdown of the (high-energy) underlying symmetry, the 
Little Higgs lagrangian below the scale Anf 7] can be written as a non-linear sigma model based on a coset 
SU(5)/SO(5) symmetry: 

£ S = I^Tr|^E| 2 , (1) 

where the subgroup [SU(2) x £/(l)] 2 of SU(b) is promoted to a local gauge symmetry. The covariant derivative 
is defined as 

2 

P M E = d„Z-iJ2 (siO^-E + YW?) + g'^BjE + EBj)) . (2) 
To exhibit the interactions, one can expand E in powers of 1// around its vacuum expectation value Eq 
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where h is the doublet that will remain light and is a triplet under the unbroken SU(2). The non-zero 
vacuum expectation value of the field (E) = Eo leads to the breaking of the global SU(5) symmetry to SO (5) 
and also breaks the local gauge symmetry [SU(2) x ?7(1)] 2 into its diagonal subgroup, which is identified with 
the standard model SUl(2) x UyiX) symmetry group. Following the notation of Han et al. |5|], we will denote 
the usual standard model gauge bosons mass eigenstates as Wjf, Zl and Al, where the subscript L denotes 
light in order to distinguish from the heavy states with mass of order /, denoted by W^, Zh and Ah- 

The standard model fermions acquire their masses via the usual Yukawa interactions. However, in order to 
cancel the top quark quadratic contribution to the Higgs self-energy, a new-vector like color triplet fermion pair, 
t and t' c , with quantum numbers (3, 1)^ and (3, l)-y 4 must be introduced. Since they are vector-like, they are 
allowed to have a bare mass term which is chosen such as to cancel the quadratic divergence above scale /. 

The coupling of the standard model top quark to the pseudo-Goldstone bosons and the heavy colored fermions 
in the littlest Higgs model is chosen to be 

£y = -Xifeijke xy XiSjx^kyu'i + A 2 /tt /c + h.c, (4) 

where Xi — {°3 > *3j *) and £yfc and e xy are antisymmetric tensors. The new model parameters Ai , A2 are supposed 
to be of the order 
of unity. 
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The linearized part of Eq. ^ describing the third generation mass terms and interactions with the neutral 
Higgs field, denoted by h°, (before spontaneous symmetry breaking, SSB), is given by: 

A = A 2 /tt' c - iXx^n^u'i + Xiftu'i - jrt h ° h0 * u 3 + h - c - ( 5 ) 

After SSB, we write h° = l/y/2(v + H), and follow Perelstein et al. [j| in defining left handed fields t$L = *3, 
tL = t and right handed fields u' 3R = u 3 , t' R = t' c to obtain 

r ( -> P \ ( —iX\v Ai/(1 — v 2 /2/ 2 ) \ / t 3L \ .,„-,. ( ,, 
C t - { u 3R t R )\ Q )\t L )~ iX ^ Hu 3rMl (6) 

-\xjHu 3R t L - jjH 2 u 3R i L + h.c. 



In order to leave the fermion mass term in its standard form we make the field redefinitions t 3 ^ — ► ~it 3 ^ and 
II — ► — ti, in the left-handed fields resulting in the following lagrangian: 

{ p \{ Ai« Ai/(1 - v 2 /2p) \(H L \ , 



F 2 TlLj.r - 

f iK 2/ 



+\i- r Hu' 3R i L + ^H 2 u' 3R i L + h.c 



Diagonalizing the mass matrix 



we obtain the usual result for the eigenvalues corresponding to the top quark t and the heavy top T which are, 
up to order 0(v/f): 

mt = TWTri v; mT = \/ A i + A 2 f- ^ 

V Ai + A 2 

In our numerical code, we will use as input the values of m-r and /, from which one obtains the required values 
of Ai and A 2 : 



From equation (|10|) . one clearly sees that there is a condition relating top masses and v.e.v.'s that these models 
impose: 

m T > 2^ c V2f (11) 

and which we incorporate in our analysis. The relevant couplings between Higgs and top quarks are obtained 
in a straightforward manner after diagonalization of Eq. and are given by (the corresponding vertices are 
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obtained via multiplication by (— i)): 



gmt = 



Ai A2 
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The relevant Feynman diagrams for the Higgs self-energy are shown in figure |T]l. 
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FIG. 1: One-loop corrections to the Higgs mass, to order w//: (a) standard top quark loop, (b) mixture of standard and 
extra top quark loop, (c) extra top quark loop with a 4-particle vertex, and (d) tadpoles with standard and with extra 
top quark loops. There are other diagrams but they are suppressed by factors of order (v/f) 2 or higher. 

The cancellation of tadpole diagrams requires that 

9Httm t + gHTTm T = (13) 
whereas the cancellation of higgs self-energy quadratic divergences implies 

9mt + 9htt + 9HT R t L + 9m R T L + 9hhu m t + 9hhtt m T = (14) 

These conditions are satisfied up to terms of order 0(v/f) by the masses and couplings listed above. 
In the simplest LH models, strict bounds on the parameters are obtained. In particular, electroweak precision 
constraints require / > 3.5 TeV 4]. However, in a recent variation on the littlest Higgs model, where a so-called 
T-parity that interchanges the two subgroups [SU(2) x U(l)]i and [SU(2) x f(l)]2 of SU(5) is introduced, can 
significantly lower this bound to / > 500 GeV ^(J . This is an important point for the phenomenology of these 
models, since a lower / implies larger deviations from the SM. Since the T-odd states do not participate in the 
cancellation of quadratic divergences, our calculation is valid in models with T-parity as well. T-parity also 
forbids the generation of a vacuum expectation value for the triplet scalar field (i.e., v' = in the notation of 
T. Han et ai.0), which is one of the causes for easing the electroweak constraints. 
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III. AMPLITUDES FOR DOUBLE HIGGS PRODUCTION 

We now turn to Higgs boson production at the LHC in the LH model, which involves the very same couplings 
responsible for the cancellation of quadratic divergences. 

Gluon-gluon fusion is the dominant mechanism for SM Higgs boson pair production at the LHC The 
amplitude for gg — > HH is dominated by top quark loops, in the form of triangle and box diagrams. Figs. |2] 
and|3]show the case for a LH model. The SM case is similar, except that Fig. |2a and all extra heavy-top loops 
are absent. We also would like to point out that T-parity forbids a term like hh<p in the radiatively generated 
Coleman- Weinberg potential. Therefore, there is no contribution of the heavy scalar in Fig. [2b and the trilinear 
Higgs coupling is the same as in the SM. 
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FIG. 2: Triangle contributions to Higgs boson pair production at LHC in a Little Higgs model. 
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FIG. 3: Box contributions to Higgs boson pair production at LHC in a Little Higgs model. 

Let us now write the expressions for the amplitude. In what follows, the external momenta p a , pb, p c , Pd are 
defined as incoming. The contribution from triangle diagrams is given by: 



iM A (g a gb -* H c H d ) = -^tSab ( gmtl(m t ) 



9HHH 



47T 3 



- gHTTlimr) 



s - M 2 H - iM H T H 



9hhh 



s-Mfj- iM H T H 



gHHttI{m t ) + gHHTTl{rriT 



where the integral /(mg) is: 



w \ _ [ H 4 n Tr U + m Q )r(4 +j>a + m Q )Y(4 + j>a + j>b + m Q )\ 

nm Q )-Jdq [q 2_ m 2 Q][{q+Va? _ m 2 Q][{q+pa+pb? _ m 2 Q] ^PaK{Pb) 

This integral reduces to the following result: 

J(m Q ) = i 4?r 2 m Q [l + (2m 2 Q - J/2)C (0, 0, s, m 2 Q , m 2 Q , m 2 Q )]e{p a ) ■ e(p b ), 



(15) 



(16) 



(17) 



G 



where Cq is the scalar Passarino-Veltman integral 12j defined as: 

A 1 

iir 2 [q 2 - m Q ] [(q + p a ) 2 - m 2 Q ] [(q + p a + p b ) 2 - m 2 Q 
The contribution from box diagrams can be written as 



iM a {g a g b ^ H c H d ) = --^Sab {g 2 Hu h{m t ) + g 2 HTT h{m T ) + h{mt,m T )) (19) 

Because the Higgs vertex is not diagonal in the (i, T) flavor, there are two types of boxes. The function I\ 
comes from boxes with either only standard top quarks or only extra heavy-top quarks in them, whereas the 
function I2 comes from boxes with both tops and extra heavy-tops. 

There are two basic box diagrams, planar (Fig. |3]a) and non-planar (Fig[S|b), according to whether the 
Higgses are adjacent in the loop or not. A planar box with a single type of quark is given by: 

Tr [(4 + mg^d + j> a + m Q )Y{j + j>a + A> + m Q ){j + j> a + j) b + j> c + m Q )} 

2 1 \( n _L n \2 — -m 2 1 \( n -X- r, -J- r>, \2 — m 2 1 \l n -X- m J- 0-1, J- n "\ 2 _ m 2 



7f (m Q ) = / d\ 



[q 2 - m 2 Q ] [(g + p a ) 2 - m Q ] [(q + p a + Pb) 2 ~ m Q ] [{q +p a +Pb+ Pc) 



m, 



xe^,(p a )e v (j)b) + \p a <-> Pb} + \pc <-> Pd} + |p a <-> Pb, Pc Pd }, (20) 
while a non-planar box is: 

4 Tr [(4 + m Q )Y l H + i> a + m Q ) (rf + i>a +lt>c + m Q )Y H + + ib + fa + m Q )] 



h NP (m Q ) = J A 



[q 2 - m 2 Q ] [(q + p a ) 2 - m 2 Q ] [(q + p a + p c ) 2 - m 2 Q ] [(q + p a + p b + p c ) 2 - m 2 Q } 

xe M (p a )ei,(pb) +{p a ^Pb} (21) 

The total contribution for boxes with a single type of quark is then 

h{m) = I p {m) +I NP {m) (22) 

We also have to compute the contribution of box diagrams with both t and T running in the loop. There are 
also planar and a non-planar contributions in this case. For the planar contribution we have: 

t p, , f ,4 Tr U + m t)r(4 + fa + rn t )Y{j + i> a +j> b + m t ) {g H T R t L ^ + 9m R T L ^)] 

I 2[ m t , mT ) J aq [q 2_ m 2 ][{q+pa) 2_ m2t][{q + pa+pbf _ m 2 ] 

x [(4+fa+fa+fa + TOT) (9HT RtL + ggt fl T^ ^)] 

[(g + p a + Pb + Pc) 2 - ra 2 T \ 
xe M (Pa)e I '(Pfc) + |p a <-> Pb} + |p c <-> Pd} + |pa <-> Pb, Pc ^> Pd} (23) 

and the non-planar contribution is written as: 

Tr [{i + m t )^{i + i> a + m t ) (g H T R t L ^ + gm R T L ^) d + fa + i>b + m t )Y] 



I? p (m t ,m T ) = J d 4 q 



[q 2 - to 2 ] [(q + p a ) 2 ~ m 2 } [(q + p a + Pb) 2 ~ mf] 
[d + j>a + j>b + j>c + m T ) (gHT R t L ^2 + gm R T L 
[(q + Pa+Pb+ Pc) 2 - m 2 T \ 

X^(PaV(Pb) + [pa ^Pb } (24) 

Accordingly, the total contribution for boxes with both t and T is given by: 

h = l[(m t ,m T ) + ^(mr^t) + I? P (m t ,m T ) + I? P (m T ,m t ). (25) 

We can then express these integrals in terms of Passarino-Veltman functions, in a way analogous to Eqs. (|16|l 
and <|17|l . We computed these transformations using the package FeynCalc 01 • Since this procedure is straight- 
forward and the final result is rather long, we did not include the expressions here. 
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IV. CROSS SECTION RESULTS 



Using the total scattering amplitude M(g a gb — > H c Hd) — M A + M a we can build the partonic differential 
cross section: 



da 1 



dfl 128tt 2 s 



yJl-4M%/s \M\ 2 . (26) 



We must point out that we have included a factor of 1/2 due to the identical particles in the final state. 
Consequently, to obtain the total cross section one must integrate Eq. 112611 over the whole 4tt solid angle. Here 
\A4\ 2 is the squared matrix element averaged over initial color and helicity states: 

1 



\M\ 2 = ^ £ IMil 2 (27) 

i,i=l,2 

where the sum is over the two physical gluon polarizations. 

We performed the calculation in the center-of-momentum frame of the gluons. In that case, the transversality 
condition of the gluon polarization vectors also implies p a ■ e(pb) — Pb ■ e iPa) — 0. Therefore we use the following 
parametrization (recalling that all 4-momenta were defined as incoming): 

Pa = (VI/2,0,0,VI/2) (28) 

Pb = (VI/2,0,0,-VI/2) 

p c = (— vI/2,0, — psin#, — pcosfl) 

Pd — (— V^/2, 0,psin#,pcos#) 

e« = (0,1,0,0) 

e (2) = (0,0,1,0) 



where the Higgs boson center-of-mass momentum is p = a/s/4 — M% . We numerically integrate the Passarino- 
Veltman functions using the package LoopTools [l4| . Finally, we obtain the pp — > HH cross section at the LHC 
by convoluting the partonic cross section with the gluon distribution function: 

a(pp^HH) = ^- J d Xl dx 2 [g 1 (x 1 ,Q 2 )g 2 (x 2 ,Q 2 ) + g 2 (x 1 ,Q 2 )g 1 {x 2 ,Q 2 )]a(gg ^ HH)6{x lX2 s - AM 2 H ), (29) 

where we used the Set 3 of CTEQ6 leading gluon distribution function with momentum scale Q 2 = s [15J . A 
K = 2 factor was included to take into account QCD corrections [lfij . 

In Fig.0|we plot the cross section for the double Higgs production process at the LHC for fixed Mt = 4 TeV, 
a Higgs boson mass in the range 150-300 GeV and for / = 500, 1000 and 2000 GeV. As expected, we find that 
the largest deviations from the SM result occurs for small Higgs boson mass and small decay constant /. In 
this sense it is important to consider models with T-p arity, where / is not required to be too large. Our results 
are otherwise consistent with the authors of Ref. where values around / = 3.5 TeV are used. 

In order to explore the dependence on the heavy top quark mass we plot in Fig. [S] we plot the cross section 
for the double Higgs production process at the LHC for fixed Mh = 200 GeV and / = 1000 as a function of 
Mt- We can see that the result grows with mr, reaching a constant limit for mr above ~ 2.5 TeV. 

We also include an analysis of the significance of the signals, which defined as: 

. £-vlh — £--&sm , orA 
significance = ^=^= , (30) 

V X • a sm 

where C is the integrated luminosity of LHC and ulh, asM are the two-Higgs production cross sections for pp 
collisions in the Little Higgs and Standard models, respectively. Assuming a luminosity of C = 10 fb _1 for the 
LHC, we obtain the results shown in Fig. EI Of course the significance scales as sfC, and the figure can be easily 
used for larger LHC integrated luminosities. 

A word of caution is necessary at this point. It is beyond the scope of this work to include an analysis of 
the simulation of detector sensitivity and possible backgrounds needed for a realistic study of double Higgs 
production and detection in the LH model, which was done in the SM in 01 • 0m study is just a first step in 
gauging how different from the SM the LH signal for double Higgs production is. 
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FIG. 4: Cross section for double Higgs production at the LHC for Mt = 4 TeV and / = 500 GeV (dashed line), 1000 
GeV (short dashed line) and 2000 GeV (dotted line). In solid line is shown the SM result. 
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FIG. 5: Cross section for double Higgs production at the LHC for Mh = 200 GeV and / = 1000 GeV, as a function of 
the heavy top quark mass Mt- 

V. CONCLUSIONS 

The double Higgs production process probes the nature of the electroweak symmetry breaking mechanism. 
This process is intimately tied to the cancellation of quadratic divergences in Little Higgs models. Here we 
have studied the reach of the LHC to probe the LH models in this way. We found that only for relatively small 
values of the energy scale /, of the order of 500 to 1000 GeV, it is possible to distinguish meaningfully the LH 
from the SM. These low values are attainable without violating the electroweak precision limits only in models 
where an extra T parity is incorporated O n the other hand, these results only mildly dependent on the 
heavy top quark mass my; while the situation is more promising for larger values of mx, it becomes practically 
independent of it for mj- above 2.5 TeV. 



Ql , ' ' ■ ■ ■ • ■ ~ ^ 
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FIG. 6: Significance of the two-Higgs signal in the Little Higgs model with respect to the Standard model, at the LHC, 
assuming an integrated luminosity of 10 fb~ . Solid, dashed and dotted lines correspond to / = 500, 1000 and 2000 
GeV, respectively. 
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